May 2001 



Solutions of the atmospheric, solar and LSND neutrino anomalies 
from TeV scale quark-lepton unification 

T. L. Yoon and R. Foot 



School of Physics 
Research Centre for High Energy Physics 
The University of Melbourne 
Victoria 3010 Australia 



Abstract 

There is a unique SU{&) ® SU{2) L ® SU{2) R gauge model which allows quarks and 
leptons to be unified at the TeV scale. It is already known that the neutrino masses arise 
radiatively in the model and are naturally light. We study the atmospheric, solar and LSND 
neutrino anomalies within the framework of this model. 



I. INTRODUCTION 



The similarity of the quarks and leptons may be due to some type of symmetry between 
them. Various theories have been proposed including the Pati-Salam theory |lj where the 
leptons take the fourth colour within a S'L r (4) ® SU{2) L ® SU{2) R gauge modelQ While 
definitely a good idea, there are some slightly unpleasant aspects of the Pati-Salam model. In 
particular, one of the main drawbacks of the Pati-Salam theory is that almost all of its unique 
predictions for new physics cannot be tested because the experimental constraints on the 
symmetry breaking scale imply that it is out of reach of current and proposed experiments. 
The problem is twofold. First there are stringent constraints coming from rare meson decays. 
These imply a lower limit on the symmetry breaking scale of about 20 TeV for the 
symmetry breaking scale which means that the heavy gauge bosons are too heavy to be 
found in even the Large Hadron Collider (LHC). The Pati-Salam model at the relatively 
low scale of 20 TeV also has great problems in explaining the light neutrino masses. The 
see-saw mechanism adopted in such models cannot suppress the neutrino mass sufficiently 
unless the symmetry breaking scale is very high. The required light neutrino masses suggest 
that the symmetry breaking scale is at least about 50 PeV ||] (1 PeV = 1000 TeV). All 
is not lost however. There appears to be a unique alternative SU{4) ® SU{2)l ® SU(2)r 
gauge model (which we call the alternative 422 model) which preserves the elegance and 
simplicity of the original Pati-Salam theory while allowing for a low symmetry breaking scale 
of about 1 TeV |J. This not only allows the unique predictions of the theory to be tested 
in collider experiments, but the theory also avoids the dreaded gauge hierarchy problem by 
not introducing any hierarchy to begin with. The theory also has characteristic predictions 
for rare B, K decays, baryon number violation as well as non-zero neutrino masses, all of 
which are naturally within current bounds, despite the low symmetry breaking scale of a 
TeV. Thus in one act of prestidigitation all of the problems aflicting the original Pati-Salam 
model are cured. 

Over the last few years, significant evidence for neutrino masses has emerged from the 
neutrino physics anomalies: the atmospheric |J, solar 0] and LSND § neutrino experi- 
ments. It is therefore an interesting question as to whether the alternative 422 model can 
accommodate these experiments. It has already been shown || that the masses for the 
neutrinos in this model typically span the necessary range to possibly account for these 
experiments. In fact viewed simply as a gauge model for neutrino masses, the theory is 
quite interesting because it provides a nice explanation for the small masses of the neutri- 
nos without any need for (untestable) hypothesis about high energy scales (which arise in 
most popular theories of neutrino masses). We will show that the theory in its minimal 
form can accommodate the atmospheric and LSND neutrino anomalies but not all three 
(including solar) simultaneously. Thus, the theory is a candidate for the physics responsible 
for the neutrino physics anomalies. However, since it cannot explain all three of the neu- 
trino anomalies, it obviously follows that if all three anomalies are confirmed in forthcoming 



* Other possibilities include models with a discrete quark-lepton symmetry which features a 
spontaneously broken SU(3) colour group for leptons. pi 
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experiments, then this would require physics beyond this model for an explanation. One ele- 
gant possibility is the mirror symmetrized extension which can provide a simple explanation 
of the neutrino physics anomalies, as we will show. 

The outline of this paper is as follows: In section II we briefly comment on the exper- 
imental situation and the various oscillation solutions to the solar, atmospheric neutrino 
anomalies and the LSND measurements. In section III we revise the essentials of the al- 
ternative 422 model. In particular we will derive the mass matrices of the fermions after 
spontaneous symmetry breaking (SSB). We will also define a basis for the fermions in which 
their weak eigenstates are related to their respective mass eigenstates via CKM-type unitary 
matrices. In section IV we explain the mechanisms that give rise to the masses of the neu- 
trinos. This includes the tree level mixing that generates right-handed neutrino Majorana 
masses f\ z>rM#(z/r) c and two radiative mechanisms, via gauge and scalar interactions, that 
give rise to Dirac masses of the ordinary neutrinos. In section V we consider the special case 
of decoupled generations and examine the possibility to obtain near maximal ul {vr) c 
oscillations within the model (which turns out to be negative). In section VI, we mirror 
symmetrise the alternative 422 model to obtain a TeV scale solution scheme for all three 
of the neutrino anomalies. In section VII we will identify how the minimal alternative 422 
model could provide simultaneous solutions to (near maximal oscillations) atmospheric neu- 
trino anomaly and the LSND measurements in the case where gauge interactions dominates 
the radiative neutrino mass generation. We conclude in section VIII. 



II. OSCILLATION SOLUTIONS TO THE ATMOSPHERIC AND SOLAR 
NEUTRINO ANOMALIES AND LSND MEASUREMENTS 



The experimental situation regarding the oscillation solutions of the atmospheric and 
solar neutrino problems has made much progress over the last few years. In this paper we 
mainly focus on the simplest possible solutions which involve maximal (or near maximal) 
two-flavour oscillations. 



A. Atmospheric neutrino anomaly 



In the case of atmospheric neutrino anomaly there is compelling evidence that about half 
of the up-going flux disappears 0. The simplest oscillation solutions which can explain 
the data are maximal z/„ — > v T or maximal Vn — > ^sterile oscillations [|10||. Despite impressive 



efforts by Super-Kamiokande the experimental data cannot yet distinguish between 
these two possibilities ||12|| . Unfortunately, this situation probably cannot be clarified until 
long-baseline experiments provide or fail to provide r events approximately in the year 2007. 



^Throughout this paper the tilde on the fermion fields is to signify that they are flavour eigenstates 
as opposed to mass eigenstates (without a tilde) in this model. See Section III for more detail. 
Note that the tilde has nothing to do with supersymmetry other than a mere notation coincidence. 
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At present the parameter range that is consistent with the atmospheric neutrino anomaly 
is roughly sin 2 26 ~ 0.85 and 



10- 3 & 5m 2 at JeV 2 ~ 10~ 2 . (1) 



B. Solar neutrino anomaly 

In the case of solar neutrino anomaly there is very strong evidence that about half 
of the v e flux from the sun has gone missing when compared to theoretical models. The 
simplest explanation of this is in terms of maximal v e oscillations, with the main suspects 
being maximal v e —>■ u a [v a is some linear combination of or v T ) JTBHIH] or maximal 
v e ^sterile oscillations |15]-[nj [|. This maximum oscillation solution to the solar neutrino 
problem can explain the 50 % flux reduction for a large parameter range: 

3 x 10" 10 eV 2 ~ 5m 2 ~ 10~ 3 eV 2 . (2) 

The upper bound arises from the lack of u e disappearance in the CHOOZ and Palo Verde 



experiments [ 19| , while the lower bound comes from a lack of any distortion in the measured 
Super-Kamiokande recoil energy spectrum [pO| , pTf , which should make its appearance for 
~ 3 x 10 10 eV 2 (traditional 'just so' region). Note that, the maximum oscillation 
solution was the only oscillation solution to predict the approximate energy independent 
spectrum obtained by Super-Kamiokande. (SMA MSW, LMA MSW and 'Just So' all pre- 
dicted some distortion that should have been seen.) The Super-Kamiokande collaboration 
has also searched for a day-night effect. However no evidence for any difference in the day 
and night time event rates were found with a 3<r upper limit of | 20| 



An- d < 0.055 (3) 

where A n _d = (N — D)/(N + D) (N = night time events and D = day time events). This 
limit allows a slice of parameter space to be excluded (using the numerical results of Ref. 

2 x 10~ 7 eV 2 ~ \5m 2 solar \ ~ 10~ 5 eV 2 (sterile) 

4 x 10~ 7 eV 2 ~ \6m 2 olar \ ~2x 10~ 5 eV 2 (active) . (4) 

Thus for both the active and sterile maximal oscillation solutions the allowed 5m 2 range 
breaks up into a high 5m 2 region and a low 5m 2 region: 

2y x 10~ 5 ~ W/eV 2 ~ 10~ 3 (high 5m 2 region) 

3 x 10~ 10 ~ W/eV 2 ~ 4y x 10~ 7 (low 5m 2 region), (5) 



^Note that maximal (or near maximal) v e — ► v sterile oscillations (and/or maximal — > v sterile 
oscillations) are consistent with standard big-bang nucleosynthesis (BBN) for a large range of 
parameters 
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where y = 0.5 for sterile case and y = 1 for the active case. These oscillation solutions will 
be tested in the near future by SNO, Borexino and KamLAND experiments. In fact the 
SNO experiment has recently announced their first results |22|] which is a measurement of 
the charged current event rate. This result, when combined with the elastic scattering rate 
obtained at super-Kamiokande disfavours the v e — > v ste riie oscillation solution at about the 3 
sigma level p2 . However both measurements are dominated by systematics which suggests 
that this result is not yet convincing but is nevertheless an interesting hint. Things should 
become much clearer when SNO measures the neutral current event rate which should allow 
v e v a t° be distinguished from v e —>■ v steT n e at more than 7 sigma. Meanwhile, Borexino 
p3| can test the low 5m 2 region by searching for a day-night effect and also seasonal effects 
fig] while KamLAND |24] can test the high Sm 2 region by searching for v e disappearance. 
Finally part of the high 5m 2 region, 10~ 4 ~ Sm 2 /eV 2 ~ 10~ 3 impacts on the atmospheric 
electron-like events and is currently disfavoured p5|. 



C. LSND data 

There is strong and interesting evidence for — + v e oscillations coming from the LSND 
experiment which suggests the parameter region 

0.2 - Sm^/eV 2 - 3 (6) 

with a rather small mixing angles sin 2 29 ~ 3 x 10~ 2 — 10~ 3 . This result will be checked by 
BooNE |IB| soon. If the LSND signal is verified then one must invoke additional (sterile) 
neutrino(s) to simultaneously explain all the solar, atmospheric and LSND data due to the 
large Sm 2 ^^ ~ eV 2 gap. Even if LSND is not verified, effectively sterile neutrino may 
still be responsible for the solar and/or atmospheric neutrino anomalies. The origin of 
the atmospheric and solar neutrino anomalies is something that only careful experimental 
studies can establish. 



III. THE MODEL 

In this paper we shall study the physics of the neutrino masses in the alternative 422 
model . Before doing so it is instructive to revise the essentials of this model (with some 
refinement). The gauge symmetry of the alternative 422 model is 

SU(i) ® SU(2) L ® SU(2) R . (7) 

Under this gauge symmetry the fermions of each generation transform in the anomaly free 
representations: 

Q L ~ (4, 2, 1), Q R ~ (4, 1, 2), f L ~ (1, 2, 2). (8) 

The minimal choice of scalar multiplets which can both break the gauge symmetry correctly 
and give all of the charged fermions mass is 
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(4,2,1), Xfl ~(4,l,2), 0~ (1,2,2). 



(9) 



Observe that the required scalar multiplets have the same gauge representation as those of 
the fermions which gives some degree of elegance to the scalar sector (although there are 
three generations of fermions and only one generation of scalars). These scalars couple to 
the fermions as follows: 



C = AiTr 



QlUlY^Xr 



A 2 Tr 



+ A 3 Tr 



Ql^Qr 



+ A 4 Tr 



Ql<P c t 2 Qr 



+ H.c, 
(10) 



where the generation index has been suppressed and C = r 2 0*r 2 . Under the SU (3) c <8>[/(1)t 
subgroup of SU(A), the 4 representation has the branching rule 4 = 3(1/3) + 1(— 1). We will 
assume that the T = — 1, I 3 r = 1/2 (I 3 l = 1/2) components of Xr{Xl) gain non-zero Vacuum 
Expectation Values (VEVs) as well as the I 3 l = —I 3 r = — 1/2 and I 3 l = —I 3 r = 1/2 
components of the (f>. We denote these VEVs by wr,l,ui,2 respectively. In other words, 



( XR (T = -1, I 3R = 1/2)) = w R , ( X l(T = -1, 1 3L = 1/2)) = w L , 

(0(/ 3 L = ~h R = -1/2)) = U U (0(/ 3 L = ~hR = 1/2)) = U 2 . 



(11) 



We will assume that the VEVs satisfy w R > u^ 2 ,w L so that the symmetry is broken as 
follows: 



SU(4) ® SU(2) L ® SU(2) R 
I (xr) 

SU(3) c ®SU(2) l ®U(1) y 

i ((/>), (Xl) 
SU(3) C ®U(1) Q 



(12) 



where Y = T + 2J3R is the linear combination of T and which annihilates (xr) (i-e. 
Y(Xr) = 0) an< l Q = hi + -j is the generator of the unbroken electromagnetic gauge 
symmetry. Observe that in the limit where wr ^> wl,ui,u 2 , the model reduces to the 
standard model. The VEV wr breaks the gauge symmetry to the standard model subgroup. 

To facilitate easy reference, we will use a — ±|, f3 — ±| to index the SU(2) l and SU (2) R 
component respectively, whereas 7 = {7', 4} is used to index the SU(4) components, where 
7' = (y,g,b) is the usual colour index for SU(3) C , and 7 = 4 the forth colour. With this 
index scheme the fermion multiplets are written as (with the generation index suppressed) 




Qr 1 





(13) 



The rationale to label some of the fermion fields in the multiplets with a tilde will be 
addressed shortly in the next paragraph. In the above matrices the first row of Ql and 
/l (Qr) is the I 3 l (I^r) = a (/3) = 1/2 component while the second row is the I 3 l (I3R) — 
a (P) = — 1/2 component. The two columns of Q L , Q R are the 7 = 7' and 7 = 4 components 
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of <S77(4), and the columns of fi are the I% R — (3 = ±1/2 components. Each field in the 
multiplets Eq. fllTf ) represents 3x1 column vector of three generations. 

As in the case of the Standard Model, the flavour states of the fermion fields are, in 
general, not aligned with the corresponding mass eigenstates. This shall necessitate the 
introduction of some theoretically arbitrary CKM-type unitary matrices into the theory. 
Without loss of generality we can choose a basis such that Dl (the left-handed down- type 
quarks) in Ql, l R (the right-handed charged leptons fields) in Q R and li (left-handed charged 
lepton fields) in fi are (almost |) mass eigenstate fields that couple with their respective 
diagonal mass matrix. The rest of the fields ip (fields with a tilde) are flavour states. 

It is instructive to list and label the scalar fields index explicitly as follows: 



4,1 \ / 4 i \ / ,1.1 ,11 



y - 



)2>2 (7)2' 2 



The electric charges of the components of Xl,r, 4> can De rea d off from Eq. (|13D because the 
scalars and fermions have the same gauge representation. 

The mass matrices (after spontaneous symmetry breaking) of the fermions are derived 
from the Yukawa Lagrangian of Eq. ( |10|) as follows: 

C (after SSB) = C E + C l + C u + C d + H.c, (15) 

where 



and 



- iC E = ElM E E R -E° L M E E° R , -iC l = (E R ) c Mii) R + \M l l R , 

-iC u = U~ L M U U R + ElM u v R , -iC d = -W L M d D R - Ej;M d l R , (16) 

M E = w R Xi, M u = (\ 3 u 2 - A 4 ui), M d = (A 4 u 2 - A3W1), (17) 



are 3x3 generally non-diagonal mass matrices for the exotic E leptons, up-type quarks U 
and down-type quarks D respectively. 



Mi = w L \ 2 



fm e 

E m M 1 (18) 
\ m T 



is the diagonal mass matrix for the charged leptons. The mass matrices Me, M u , M d can 
be diagonalised in the usual way (biunitary diagonalisation) to obtain the diagonal mass 
matrices, e.g. Me = U^MeV etc. so that the flavour states ip are related to their mass 
eigenstates ip via unitary matrices: 



§ Strictly, II, l R are approximately but not exactly mass eigenstates because of the mass mixing 
between l R and E~[ [see the forthcoming Eq. (p])]. This means that the true charged lepton mass 
eigenstate fields have the form Z|™ e = l L + 0{^M±)E L , l% ue = l R + 0(^)E R . 
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E R = V*E R , E L = U^E L , 

u r = y r u r , u l = y}u l , 

D R = K'D R , (19) 

and Dl = -Dl, II = Il,Ir = Ir with our choice of basis. Note that approximately the same 
unitary matrix U (V) relates both E®, E^ (E R , E R ) to their weak eigenstates E®, E£ (E R , 
E R ) because the mass matrices are approximately SU{2) L (SU(2) R ) invariant. 

The gauge fields in the 15 representation of SU(4) has the branching rule 15 = 8(0) + 
3(— |) +3(|) + 1(0) under SU(3) L (g)U(l) T . We identify the colour octet 8(0) as the gluons of 
the usual SU(3) C colour group, 3(— |), 3(|) colour triplet lepto-quark gauge bosons W, W* 
that couple the 7 = 7' components to 7 = 4 component in Xl,r, Ql,r- Note that there is 
also a neutral gauge boson B' corresponds to the singlet 1(0). The matrix = is the 
usual CKM matrix (as in the Standard Model), whereas Y R K'^ = Kr is the analogue of the 
CKM matrix for the right-handed charged quarks in SU(2)r sector. The matrix K' is the 
analogue of the CKM-type matrix in the SU(4) sector pertaining to lepto-quark interactions 
mediated by W, W* . It was shown in Ref. ||[| that the main experimental constraints on 
this 422 model come from rare B and K decays such as K° — > ^e^^B — ► /i ± r :F etc. 
depending on the form of K' . Remarkably, the symmetry breaking scale could be as low 
as a TeV without being in conflict with any experimental measurements. Indeed we will 
assume that the symmetry breaking scale is in the interesting low range: 

0.5(1.0) TeV ~ M Wr (M w ,) ~ 10 TeV, 45 GeV ~ M E ~ 10 TeV, (20) 

which is well motivated because it avoids the gauge hierarchy problem and it also allows the 
model to be testable at existing and future colliders (such as LHC). (Note that the lower 
limit on the mass of the E lepton arises from LEP measurements of the Z° width). 

Apart from laboratory experimental bounds, there are also astrophysical ones on My/ R - 
In particular Ref. [p7j-p9| have argued that based on 'energy-loss' argument, Supernova 
1987A excludes a range of values for Mw R and the Wl — Wr mixing parameter, ( Q (for 
m UR ~ 10 MeV): 



C< 10 



-5 



(0.3 - 0.5) ~ - 22 -40 in the limit C -> 0, (21) 



which seems to marginally rule out the range of Mw R as assumed in Eq. (P0[). However, 
while there will undoubtly be effects for supernova coming from the additional gauge bosons 
in our model, we should nonetheless keep in mind that the modeling of core collapse of 
supernova is generally plagued by theoretical, observational as well numerical uncertain- 
ties. For example, Berezinsky [R0| has recently emphasised that the reasonably successful 



"In our case here, £ ~ — ~4xl0 5 — 1 x 10 4 , where = gL9RUiU2, see forth coming Eq. 
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description of SN1987A is somewhat surprising given that it was assumed that the presu- 
pernova protostar was a non-rotating red super giant, while it appears that it was actually 
a rotating blue super giant. Also, Turner [STlj] pointed out that there exists uncertainty in 



the theoretical model for the hot core of a core collapsed supernova, which itself depends 
critically upon the equation of state at supernuclear densities (which is a state-of-art prob- 
lem in nuclear physics). The criterion of neutrino luminosity from SN 1987A, Q a ~ 10 53 
erg s _1 , a key ingredient in obtaining the SN 1987A bound of Eq. ([H]), according to Ref. 



3j| is also subjected to question. Raffelt and Seckel [^7]] pointed out that SN 1987A bound 



on the right-handed and other light exotic particle interactions could be uncertain up to as 
much as 2 order of magnitude, in which case could render the apparent 'conflict' of the range 
assumed by Eq. ( p0|) and the SN 1987A bound of Eq. (pT|) be alleviated. While keeping 
in mind the possible astrophysical implications of a low symmetry breaking scale Myy R ~ 1 
TeV, we now continue with our exploration of the possible phenomenology of the alternative 
422 model. 

The SU(2) LR charged gauge bosons Wf; R couple to the fermions via the interaction 
Lagrangian density 

lC „ = J^p£ W + KLDL+ ^ W +i L + El W +Ei + {E^WUE R ) C + H.c] 

+ % Wr W£k r d r + w R ^ + i R + TW c v^ Wr- k + Je^y v* WPl 

+ H.c.]. (22) 

The charged gauge interactions are of our interest because they will give rise to radiative 
Dirac mass terms to the neutrinos in this model as we will now discuss. 



IV. NEUTRINO MASS 

With the model as defined in Section III, the ordinary neutrinos vl are massless at 
tree level because the v L states do not couple to any VEV [see Eq. (UDO- However, 
the neutrino masses are non-zero in the model because there are 1-loop (and higher order) 
Feynman diagrams which contribute to their masses. In other words the masses of the 
neutrinos arise radiatively in the model. In particular, as we will see later, a Dirac mass vtld 
and a &l — {El) c mass mixing term m V E will be generated as mass corrections at 1-loop level, 
meanwhile the v R states gain Majorana masses at tree-level by mixing with E leptons. We 



^ To generate tree- level neutrino mass we need to either admit the gauge invariant bare mass term 
f^baref h{f l) c into the Lagrangian density in Eq. (|io|) or add a new Higg A ~ (4,2,3) into the 
Lagrangian density via the coupling AAA^Q/j/i+H.c.. By developing a VEV, A can generate a 
Dirac mass term ULm£>D R . We argue that, since the scale of m\, are is completely independent of 
the weak scale, the assumption that mj, are <C M wea ^ is surely an interesting possibility. In view of 
this plausible assumption we have set mi, are = 0. Meanwhile, adding an additional scalar multiplet 
such as A spoils both the simplicity and elegance of the model. 
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will elaborate these mechanisms in more detail in the following subsections. In this paper 
we will be working exclusively in the t'Hooft-Feynman gauge. 



A. Tree level Majorana mass matrix Mr 

At tree level, mixing between vr with E\ R generates right-handed neutrino Majorana 
mass Mr. In the mass eigenstate basis (for the E° LR ) defined in Eq. (|T3|), the tree-level 
Lagrangian density of Eq. (p~5|) becomes 



-iC 



( v L (u R y El (E 



M 



M d 
WM d K' 



( (vl) c \ 

(E° L y 

\ (K) J 

( D R 
{ Ir 



Dl El 

+ElM E E R + l L M t l R + U L M U U R + H.c, 



(23) 



where 



M 



/0 \ 

y r m u yIu M t V 

(Y R M u Yluy -M E 

\0 (M t Vy -M E J 



[12 x 12 matrix). 



(24) 



Md, M u are the 3x3 diagonal mass matrices for the down-type quarks and up-type quarks 
respectively, whereas Me is the 3x3 diagonal mass matrix for E leptons after biunitary- 
diagonalising Me = wlXi- 



( M, 



M f 



um f v ] = u 



Ei 







M E2 
V o M Es 



(25) 



The matrices U, V describe the relation between the weak and mass eigenstates of the E 
leptons and can be determined from Ai [see Eq. (|T9D1. 

One can block diagonalise the 9x9 matrix in the lower right sector of M by a similarity 
transformation |32| using the approximately orthogonal matrix 



h P P 
~Pp h 



(9x9 matrix) 



where 



P P 



Y R M U Y[U M^ 



-Me 1 



-M 



E 



. (3x6 matrix) 



(26) 



Block diagonalisation casts M into the form 
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M 



/ \ 
M R 




V o o - M '* 



(27) 



where each '0' is a 3 x 3 matrix of zeros, and 



M R ~ (MiV M^U^YlMuYP} + (MiV Me^YlMuY^ , (3x3 matrix) 
( M E ' 







. (6x6 matrix) 



(28) 



In the limit Me ^> Mi, M u , the states vr are decoupled from the E leptons. In the special 
case of decoupled generations (e.g. Yl = Yr = U — V = I), M reduces to 



/o 











\ 








m qu 


m 












-M Ei 




\o 


mi 










M 



where i, q u and / index the generations. In this case, Eq. (|28D reproduces the result 



(29) 



M, 



f! 



M, 



(30) 



as obtained in Ref. 



B. Radiative correction to M due to gauge interactions 

At the 1-loop level the gauge interactions from the charged SU(2) L gauge bosons W L 
and SU (2)r gauge bosons Wr give rise to vl{Pl) c Majorana mass tum, vlVr Dirac mass mo 
and Pl{E < i) c mass mixing term m vE - The Dirac mass mp arises from the gauge interactions^] 

+ ^fWAfe + H.c. (31) 

which leads to the Feynman diagram: 




^The (unphysical) Goldstone boson contributions will be evaluated together with the (physical) 
scalar contributions in Part C of this section. 
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Figure 1. Dirac mass generated by gauge interactions leading to the mass term Ul^r. 

Note that since these interactions do not mediate cross generational mixing among the 
neutrinos, the Dirac mass matrix mo is strictly diagonal. As calculated in Ref. ||, for each 
generation, 



9r9l \? , (Ml R . _ 



where /i 2 = gL9RUiU 2 is the Wj, — mixing mass. Note, that the lower limit of jj 2 is 
simply zero since /i 2 vanishes in either limits of U\ — » or tt 2 — > and the theory remains 
phenomenologically consistent in this limit. It was also shown in the Ref. that 

/2 2 < 1 M^ L m b (33) 



M 2 ^ 2^3 M^ K m t ' 

which is not a strict upper limit, but rather an approximate condition to avoid fine tuning 
0. In view of this rough upper limit, it is convenient to write in terms of a parameter 
< r] < 1 defined such that 

mo = rriik g = mirjS (34) 

where 

If we take the reasonable range of 0.5 TeV < M Wr ~ 10 TeV, then S typically spans a range 
of 

10~ 6 < S ~ 10" 9 . (36) 

Thus, the gauge loop contribution to m D is proportional to rri[ and is naturally light because 
of its radiative origin. 

The mass mixing term m U E arises at 1-loop level via the gauge interactions 

^ L W; L El + ^E^W^H + H.c. (37) 
leading to the Feynman diagram: 



To avoid the equality 



§§As shown in Ref. 0, the limit of Eq. (63) comes from 4^ < mt. 
nib = m t we will need the scale of ui,u% be separated by a hierarchy, else we will need to fine-tune 
the Yukawa coupling constants X^. 
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Figure 2. 

mass term m vE . 



vl(E l ) c neutrino mixing term generated by gauge interactions leading to the 



Note that in contrast tome, the involvement of the matrix V in the interactions Eq. ( 37|) 
may mediate cross generational mixing. In the special case of decoupled generations, m vE 
has been calculated earlier by Ref. |J as 



m vE = M E 



9r9l / fir 



r]M E S. 



m 2 Wr . 



In 



m 2 Wl , 



+ 



M 2 E - M'W L 



Mllnf^a 



M 2 E -M^ R 



(38) 



In addition to mu and m uE , this model also generates a vl(vl) c Majorana mass term rriM 
at 1-loop level by charged gauge interactions. In the case of decoupled generations, 



m M = m l m qd M J 



9r9l / H 



M w R , 



In 



M w R 



In 



m 2 Wr -mi 



M 2 — M 2 

m W L m E 



V ^4 E - b 



(39) 



is generically tiny compared to and m vE . 



C. Radiative correction to M due to scalar interactions 

In the previous subsection we have shown that at 1-loop level the gauge interactions 
give rise to Dirac mass uid, mass mixing correction m uE and also uim- Nevertheless, in 
this model this is not the only way mass corrections could arise. Other than the gauge 
mechanism discussed in the previous subsection, the Higgs sector also contains interactions 
that could generate Dirac mass for the neutrinos as a radiative mass correction. The relevant 

4-1 

interactions involve the negatively charged colour singlet scalar \i r '■ 
- iC 4 i = -iC 4 i - iC 4,-1 

x 1 Xl 2 X R 1 

4-1 4,-i 

= ^-E= L M E v\v L r - ^—v-RM^^y 

w R w l 

4-1 4-i 

1 2 '2 



ElM E V\l L ) c - ^^l R M l v L + H. c. (40) 



Wr W L 
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The first two terms of Eq. (f40|) will give rise to a Dirac mass correction mo with E~ as 
propagator (see Fig. 3), whereas the last two terms will give rise to a mass mixing term 
m U E with charged leptons as propagator. Therefore in comparison to the mass m V E can 
be safely ignored. 



E K m _ E 



E 



L 



> < X < > 



4,-1/2 - . v/ " 1 „ 4,-1/2 

XL x 

^3<(|» 



Figure 3. ^lz^j Dirac mass generated by scalar interactions. The cross in the scalar 
propagator is the perturbative mass mixing term that can be shown to vanish in either limit 
of u\ t 2 —> 

The mixing between xl — Xr is effected solely by the mixing term of the form 

V 3 = M X \$T 2 XR + M'xU C T 2 XR + H.C. 
t r 1 4,— i* , TT 

^V m =m m x L 2 Xr + H.c; 

m 2 m = -(Mm + M'u 2 ) (41) 

in the Higgs potential. 

The charged gauge bosons W LR acquire mass by eating two colourless, charged would-be 

Goldstone bosons (let's call them Gj^ and Gf;) which are linear combinations of Xlr S 0^ 
and leaving behind two physical charged Higgs (which we will call H^ 2 ). The fields 

4,-i 

Xlr appearing in Fig. 3 are linear combination of these would-be Goldstone bosons and 

physical Higgses. We shall work out the linear combinations of G LR and Hf 2 in Xlr so 
that the Dirac mass arising from these ^-interactions can be calculated. 

The Goldstone bosons associated with each spontaneously broken symmetry are given 

by 

G = $ T T a A a , (42) 

where $ are the Higgs scalars of the theory, T a the generators of the broken symmetries 
and A a the vacua. Utilising Eq. ((4^), the would-be Goldstone bosons G L and G R can be 
identified by taking T a as the charged SU(2) L R generators, r±. The two states which are 
orthogonal to G^ and G^ will be related to the physical Higgs Hf 2 with masses denoted by 
Mh 1} h 2 - We will first work in the limit u\ — > since we expect U\ <C u 2 (or u 2 C «i, which 
is analogous to u% <C u 2 in the following analysis). In this limit it is possible to show that 
(see Appendix) the weak eigenstates are related to Gf; R , H^ 2 via the unitary matrix U g as 
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4,- 

Xr 

4,- 



V0" 



" 2 ' 2 



J 



Gl 



(43) 



where the matrix U g is (in the limit u\ — > 0) 



17, 





W R 





u%/w R 




Ni 







w L 







N L 




— U2 





M2 




Ng 









«2 









\ 


u\/w L 




N 2 









"2 


J 


N 2 



fU n U 12 U 13 U 14 \ 
U21 U22 U23 U24 
U31 U32 U33 U34 

\u 41 c/ 42 u 43 U M J 



(44) 



with normalisation constants 

N 2 = w R + u 2 2 , N 2 L = w 2 L + ul 



11' 





, Nl = u\ 


f 1+ (-)l 


L \w R J . 




L V^l/ J 



(45) 



Referring to Eqs. 



we observe that 



4 -2 



Xl 



wr 



w 2 



--G-, 



R "f 
tl-2 



HI 



w\ + u 2 2 



-.m 



Thus it is clear that the x-l°°P correction vanishes in the limit u\ 
implies that 

hm V m = lim V m = 0, 



(46) 



0. Self consistency 



(47) 



which we also prove in the Appendix. In short, the x-l°°P contribution to the neutrino 
Dirac mass vanishes in the limit u\ — > (or analogously 112 — > 0). Recall that the gauge 
contribution to Dirac mass also vanishes in this limit. The expected hierarchy u\ <C U2 (or 
m 2 C Mi) thus ensures a small gauge and scalar contribution to the neutrino Dirac mass. 

We now investigate the case of small U\ 7^ (with u\ <C 112). In this case we can treat 
the U\ term as a perturbation that induces small mass mixing term Vsu\ that will couple 

4,-i 4,-i 

Xr 2 to Xl 2 j thus leading to neutrino Dirac masses. This argument also holds true if we 
interchange Ui «-> u 2 . Indexing the would-be Goldstone bosons and physical Higgs fields as 

S T a ={G~ Rl G- Ll H^H^), a = 1,2,3,4, (48) 
we can write the linear combination of G L R and H l 2 in Xlr 2 compactly as 



4,- 
Xr 



X U la >S a >, Xl 2 = X] U ^Sb- 

a'=l,3 6=2,4 



(49) 
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The explicit expression of Ui a >, U 2b can be read off directly from U g in Eq. ([44]). The mass 
correction in Fig. 3 is now a summation of all the contributions from the approximate mass 
eigenstates diagrams. Now we could evaluate the Feynman diagram in Fig. 3 by treating the 
cross insertion as a perturbation that gives rise to a vertex factor = — {Mu\ + M'u 2 ). 
In the special case of decoupled generations, the mass correction is calculated to be (for first 
generation) 



with 



where 



m D = k x m e , 



Ml A 
16tt 2 wlWr 



(50) 



(51) 



a'b 



A = Y / A a > b U la ,U 2b , 

,b 

Mu 1 + M'u 2 s 



Ml 



Ml 



M 2 S , ln(^ 



•Mi 



M 2 E ln( 



Ail 



Ml 



Ml 



Ml 



Ml 



(52) 



This quantity A can be greatly simplified in the limit w 2 R 3> u 2 ^> w\, u\. In this limit G R 
and H 2 will dominate the loop: 

'Mui + M'u 2 



A 



Mw R ~ Ml 



Ml 2 ln( 



Ml 



Mf, 



MlM^k 



M 2 H2 - M 2 Wr 



M 2 Wr - Ml 



(53) 



where we have set the mass of G R to the mass of M\y R since we are working in the t'Hooft- 
Feynman gauge. Putting in reasonable limits for the parameters (Mw R ~ 0.5 TeV, 50 GeV < 



M F 



10 TeV, 1 GeV < w L < 200 GeV), we find 



u x u 2 M\ 
w l w r M Wr 



~ k x 



< u x u 2 



W L W R 



(54) 



Thus we can now constrain the relative contribution of the gauge and x-loop diagrams. We 
find that 



M% w r < k x < w R 



M, 4 



Using jjjte- £ 1Q- 2 and 10 



< Wr < 
r^/ — ^ r*-/ 

W L 



Wr W L kg 

10 4 , we have 



w L 



(55) 



10 



-7 < k x < 



10 4 



(56) 



Thus there is a range of parameters where gauge loop dominates (k g ^> k x ) and range of 
parameters where x _ loop dominates (k g <C k x ). There is also a troublesome intermediate 
region with k g m k x which we will for the most part ignore. 
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V. u L - v R MIXING IN DECOUPLED GENERATIONS 



Having derived the various contributions to the neutral lepton mass matrix M, we now 
examine the scenario where we assume that the mixing between generations is approximately 
negligible. Referring to the tree level mass matrix M of Eq. (p4[), v R gains a Majorana mass 
from its mixing with the E leptons while is massless at tree level. Including the 1-loop 



mass corrections, the mass matrix M is 



M 



/ 



in 



D 



m D 




m uE \ 

4 



Y R M U YIU My 
m\ E (Y R M u Yluy -M E 
(M { V)+ -M E 



(57) 



The neutrinos (both left-handed and right-handed ones) will approximately decouple from 
the E leptons since the latter are much heavier (recall that Me ~ 45 GeV from Z° width). 
The effective Lagrangian density for the mass matrix of the neutrinos after decoupling from 
the E leptons is 



eff 




+ H. a, 



(58) 



where the matrix M v is as given by 



M v 



m M rri 



{m' D y M. 



D 



with 



m 



D 



mn + myEM^V^Mi. 



(59) 



(60) 



M R is given earlier in Eq. (p8|) . In the see-saw limit where the eigenvalues of M R are much 
larger than the eigenvalues of m' D , Eq. (|58"D becomes 



see— saw 



where 



C 



rriL — m M 



1 1 

-v L m L (i> L ) c + -(v R ) c M R i> R + H.c. 



(61) 



m' D M^ 



D) 



fm n 
m 2 2 
\ m 33 

/ mu 
m 2 2 
V m 33 



I Mi 
M 2 
\ M 3 



(62) 



with mu and Mi (i = 1,2,3) as the eigen masses of the mass matrices mo and M R re- 
spectively. In Eq. (|62"D we have assumed that the matrices M R and m' D are diagonalised 
with 
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Mi 
M R = U R [ M 2 | 4, m' D = Vj, 
M 3 



/mn 

m 22 | Vj. (63) 
V m 33 



Generically is always tiny in comparison to the other contributions to mj, in all gener- 
ations and shall be dropped hereafter. 

In the case of decoupled generations the mixing angle between i>l with vr can be deter- 
mined from the matrix M v : 



taa20„ = -— (64) 



M, 



r 



The eigen masses of M v are 



i(M s± ^Ti^)^( 1± -L-) (65) 
and the mass squared difference is simply 

h{ = Jl = te)'_' . (66) 
y cos 20„ \ M E J cos 2^ v ; 

Under the assumption of decoupled generations, the only way to solve atmospheric neutrino 
anomaly is via — > (v^rY oscillations. However, since experimentally we know that 
$ m atm ~ 10~ 2 e V 2 , sin 2 20 atm ~ 0.8, this implies 

M E = 2mq - mi Z 70 TeV (67) 
' Smitrr, cos 29 „ 



even with the lowest possible fermion masses m qu = m u , mi = m e . Thus it is not possible to 
solve atmospheric neutrino anomaly with decoupled generations while keeping the symmetry 
breaking scale in the interesting range ~ few TeV. The above result applies to both cases 
where m' D is dominated by the gauge loop or by the x-loop. 



VI. MIRROR SYMMETRISATION OF THE ALTERNATIVE 422 MODEL 

As concluded in the last section, decoupled generations cannot accommodate the atmo- 
spheric neutrino anomaly within our assumption of a low 422 symmetry breaking scale ~ 
few TeV. Nevertheless, we can still have near maximal active-sterile oscillations at the TeV 
scale if we mirror symmetrise the alternative 422 model in the spirit of the Exact Parity 
Model (EPM) [j3l|| . These models allows parity to be an unbroken symmetry of nature 
which it turns out offers a framework to understand the dark matter, neutrino anomalies 
and various other puzzles. (For a review of the evidence see f3~5| .) It is well known that 
in mirror matter models maximal mixing occurs naturally because of the unbroken parity 
symmetry. In the 'mirror symmetric alternative 422 model', will oscillate maximally to 
its parity partner {v' R ) c . Mirror symmetrisation can thus provide a mechanism to obtain 
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maximal mixing whereas the naturally tiny masses (hence mass squared difference) of the 
neutrinos are radiatively generated and naturally small in this TeV scale model. 

In the mirror symmetric alternative 422 model, the symmetry group is extended to 
SU{A) <g> SU{2) L <g> SU{2) R <g> SU(A)' ® SU(2)' L ® SU{2)' R . Each family has a mirror partner, 
which we denote with a prime. Considering and Ur, v' r and v' L will be included in the 
particle content enabling parity to be an unbroken symmetry. Under the parity transfor- 
mation, vl <->• "iov'foVR <-> 7o^l (as well as x —* — x, of course). The see-saw Langrangian 
density of Eq. fl6T|) becomes 



/see— saw 



1 

2 



where 



/0>l) c \ 



R 



+ H.c, 



(68) 



M' 



( 





m! 




o \ 




m! 












(m'nV 





Mr 





\ 





K) 1 








(69) 



m' is a parity invariant mass mixing term that mixes the ordinary neutrinos with the mirror 
neutrinos. In the minimal mirror matter 422 model, i.e. with only the Higgs fields (f>, xl, Xr 
and 0', x'l-i Xr> there is no mirror-ordinary mixing Yukawa coupling to generate m! . This is 
quite unlike the case of the EPM model where ordinary-mirror Yukawa coupling [y'^^v^ 
exists because vr and its parity partner v' L are gauge singlets. One could consider the 
mirror-ordinary mass mixing term as a dimension-five unrenormalisable bare mass term 
^of > R<i ) 'fL4'/Mh so that the mirror symmetric alternative 422 model becomes effectively a 
remnant resulting from the breakdown of some unknown physics at a much higher scale Mh- 
Alternatively, the required mirror-ordinary mass mixing term can be generated if an 
additional Higgs scalar p a/3,a p ~ (2, 2) (2', 2') exists. The gauge invariant Yukawa coupling 

C P = X 5 hpr 2 f' R + \Jlp c T 2 f R + H.c. (70) 



will generate a mirror-ordinary mass mixing term as p develops a VEV. Because p can couple 
to (jxft' in Higgs potential it can easily gain a VEV in the right orientation. 

To see this more transparently, let's take a look at the part of Higgs potential containing 
only the scalar field p, 



V{p) = M 2 p p^p + mp^tfp + m p {p c ) W + 0{p 3 , p 4 ) + H.c. 



(71) 



As 0, 0' develop VEVs at (0) = (0') = ( ^ ^ ) ' ^ e ^ rmnear term will induce a linear 
term in p. Because of this, minimising V(p) with respect to p^ 



dV 

— = M 2 p p + m p (0)(0') + m;(0' c )(0 c ) = 



(72) 
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will induce the VEVs 



(P) = ^ (73) 



to the components of p a P> a 'P' . The VEV for the component 



( p )(12,2'l') = ^(12,2'1') = ( m P + m p) M l M 2 ^ 

will generate the desired ordinary-mirror mass mixing term in C p : 

C p = v l wlv' r + H.c, (75) 

where 

m' = -A 5 (p)( 12 ' 2 ' 1 ') - A 6 (p c )( 12 ' 2 ' 1 ') = ^K + ^) (As + Ae) . (76) 

In any case the mirror-ordinary mass mixing term m' is still a free parameter of the theory. 
In the see-saw limit ul, (v' r ) c decouple from u' Ll {vr) c in Eq. (|68f) and the Lagrangian density 
of light neutrinos is given by 

Cl = U^ WJ C )M" + H.c, (77) 



2 s V v'r 

where 



- , -m' D M R l {m' D y w! \ , . 

M ~ 1 w! -m' D M R \m' D y J' [ ^ } 



The mass matrix M in Eq. ([Tg) describes maximal vl — > (^) c oscillations with eigen 

masses 

m± = -m' D M R l (m' D y ± m'. (79) 
In the limit of small intergenerational mixing the mass squared difference is simply 

2 _ 4m'mg _ 2m'miM Ei S 2 r) 2 

T ~ M R ~ m qu ■ m 

5m 2 ± can be identified as 5m 2 atm (second generation) and Sm 2 olar (first generation) with max- 
imal oscillation solutions, thereby explaining the atmospheric and solar neutrino anomalies. 
If small intergenerational mixing (parametrised by 9 and 0) between the first and second 
generations is included, the LSND experiment Sm 2 -^^ can be identified with 

SrrfiSND = M e+>i>+ l = \ m l e+ ~ m X + \ (81) 

and 

(sin 29 + sin 20) 2 ~ 3 x 10~ 2 - 10~ 3 . (82) 

Thus, the mirror symmetrised extension of the alternative 422 model can explain all the 
three neutrino anomalies without any physics beyond the TeV scale. The alternative 422 
model provides the neutrino masses while mirror symmetrising it provides the maximal 
mixing between each ordinary and mirror neutrino flavour. 
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VII. MIXING BETWEEN GENERATIONS: NEUTRINO MASS DOMINATED 

BY GAUGE INTERACTIONS 



In the previous two sections we have examined the case where the mixing between gen- 
erations was small. Under this assumption the minimal alternative 422 model could not 
accommodate the large mixing required to explain the atmospheric or solar neutrino prob- 
lems. However we also shown that the mirror symmetrised extension could explain all three 
neutrino anomalies. We now return to the minimal alternative 422 model and examine the 
alternative case where mixing between generations is large. From the discussion in the previ- 
ous sections, recall that there are two different ways in which neutrino masses are generated 
in the model, namely via the gauge interactions and via the (scalar) ^-interactions. In ei- 
ther case, we see that the Dirac masses are proportional to the charged lepton masses mi 
[see Eq. (0), ([SO])]. The proportional constants k x and k g are poorly constrained however. 
Our ignorance of their relative strength does not permit us to tell which mechanism is the 
dominating one. Though it remains a possibility that both contributions may be equally 
contributive, we shall not treat this general scenario to avoid complications. We will focus 
our attention only to the limiting case where the gauge interactions are assumed to dominate 
over the scalar interactions, i.e. k g ^> k x . 



A. Two generations maximal mixing and 'lop-sided' Mr 

In this subsection, we shall analyse the limiting case where the radiative contribution 
to the neutrino mass from gauge interactions dominates the contribution due to the scalar 
interactions, i.e. k g ^> k x . Our strategy is to look for a mechanism in the 2-3 sector that 
provides a near maximal — > v t l oscillation solution to atmospheric neutrino anomaly. 
We then generalise it to the case of three generations. 

Recall that the effective mass matrix for the light neutrinos is given by [see Eq. 
(|62|)1. Knowledge of mi allows us to work out the mixing angles and 5m 2 , and to make 
contact with the neutrino experiments. In general, to calculate m^ we require knowledge 
of the Yukawa coupling matrices Ai, A 4 . Our purpose here is to identify simple forms of 
Ai,...,A4 which enable the model to accommodate the neutrino data. Let us first look at 



the special case of just the 2-3 sector. As discussed in Ref. ||33|| , we can obtain maximal 
VpL —* v t l oscillations if m' D is approximately diagonal and Mr being 'lop sided', meaning 
that the off diagonal elements are much larger than the diagonal elements of the 2-3 mass 
matrix. 

The (tree level) two generations mass matrix Mr in the 2-3 sector, which we parametrise 

by 

M "={ZZ) < 83 > 

can be easily worked out from Eq. (|28|) . Note that the form of Mr2 is intimately related to 
the Yukawa couplings Ai as well as the right-handed CKM-type matrix Yr. On the other 
hand, quite independently, the radiative correction to the neutrino mass m' D can be read 
off directly from Eq. (^). Cross generational mixing of the Dirac masses is possible due to 
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the dependence of the matrix V. However, in the limiting case of diagonal V, the matrix 
m V E, and thus m' D , also becomes diagonal. In this case we can approximate m' D by 
. The Dirac masses of m' D ~ mp, namely 7^22, ^33, are simply given by Eq. (0), with 



m i = Tn^m T respectively. 

Let us parametrise the unitary matrices [as introduced in Eq. (|T^)] as 

u= /cos9 u -sm0 u \ v= /cos9 v -sin9 v \ y = /cos9 y -sm9 y \ , , 

\ sin 9 U cos 9 U ) ' \ sin 9 V cos 9 V ) ' R \ sin y cos 9 y ) ' 

The simplest way that we have found to obtain lop-sided Mr in our scheme (modulo certain 
'permutations' which we will discuss later) is to note that when Ai is approximately diagonal 
(i.e. 9 U>V small)P^|, then 

sin fly m M m e m T mt 



\\m^M R2 — ( sinjVrymf 2 m T m t 2 cos Q v m T mt ^ ) ■ (^) 



-*0 



M £2 M B3 M B3 

If we let Mr 2 be diagonalised by the unitary matrix 



sin 9r cos 9r 
we find that the mixing angle 9r behaves like 



m T m t l 'Me 3 

i hm o tan2^ = - _ ^. Q _ ■ (87) 



Mb 2 a/e 3 > wo "y 

This mean that, assuming no accidental cancellation, we will obtain |tan2#^| ^> 1 in the 
limit 

7T 

9 u ,v — > 0, 9 y — > — , (88) 

with a pair of approximately degenerate eigen masses 

/ m u m c m T mt\ , s 

M2 ^ T (^f-i^j" Tm - (89) 

Since m' D is kept diagonal due to the vanishing therefore maximal mixing in the left- 
handed neutrinos (2-3 sector) is realised. This can be seen from the effective mass matrix 
m L of Eq. (H): 



***In the case where the radiative correction is dominated by the gauge interactions m v E will 
contribute to m' D [see Eq. (60)]. In the special case of decoupled generations, m' D = + 
muEmi/MEi- Ref. || established the estimation of mo ~ "Ve^z/M^. Assuming no accidental 
cancellation, we could approximate m' D ~ mo- 

ttt Yoi simplicity we have set the CKM matrix Yl = I. 
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1 m 22 m 3 3 , . 

mL ^{io) -^r- (90) 

The effective mass matrix Eq. (|90| ) indicates that, in the presence of small intergenerational 
mixing, the left-handed active neutrinos v^i, i) t l are approximate maximal mixture of almost 
degenerate mass eigenstates. We denote these eigen masses as 

±ri 2 S 2 

™2, m 3 = mc — . (91) 

In the case where the limits of Eq. fl88|) are only approximate, the diagonal entries in 
Mr2 (and hence mi) shall in effect be replaced by some tiny values e instead. It is natural to 
conceive that e could split the eigen mass degeneracy to an order of 10" 3 eV 2 ~ |5^„ im | = 
\ m i~ m 2\ ~ ^ 2 e ^ 2 - However, if this splitting is to be natural, the size of 8m\ 2 = m\ — m\ 
must be smaller than m\,m\ themselves, i.e., 



m 



2 _2 
3' 



m\ > ID" 2 eV 2 . (92) 



Estimating m\ 3 [using the range of S from Eq. (36)] 



( A M .v. m 



Thus, we see that although near maximal — > v TL oscillations are achieved via the ansatz 
of Eq. ( [38l) the 8m? is not in natural compatibility with 5m 2 atm in this particular situation. 
This means that maximal oscillation solution to atmospheric neutrino anomaly is not accom- 
modated in this case. However we will show later that a permutation on the up-type quark 
masses in m^ could be performed to obtain compatibility with the consistency requirement 
of Eq. ( ]92|) . It is more convenient to discuss how the permutation (and its rationale) of the 
up-type quark masses could obtain a compatible range of Sm 2 between this scheme and the 
experimental values by including the first generation neutrino into the picture. We shall do 
so in the following subsection. 

Before we proceed to the three generations case, it is worth to comment on the form 
taken by under the ansatz of Eq. ([38|). Since mo is approximately diagonal under this 
ansatz, near maximal v^l — > v TL oscillations should originate from Mp® of the right-handed 
neutrino sector, which is in turn related to the form of and the mass matrix Ai. In the 
limit of this ansatz, Y^ is off-diagonal in the 2-3 sector (assuming that the first generation 
is approximately decoupled from the 2-3 sector), 

(1 0\ 

Yj ~ 1 . (94) 

If we assume a left-right similarity so that 

K R ~ K L ~ I, (95) 
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then this means that K' = Y^K R also has the form 



K' 



(I 





Recall that, as discussed in detail by Ref. 
coloured gauge bosons W'^ 




Ft. 



(96) 



the SU(4) gauge interactions involving the 



C 



9s 
V2 



DrW'^K'Ir + H. c. 

IS 



(97) 



could mediate lepto-quark transitions. In Ref. it was shown that K' must be 

non-diagonal to avoid contributions from K° — > /i ± e :p decays. In that case the pri- 
mary constraint on the SU(4) symmetry breaking scale M\y' is from rare B decays, 
Bq — > /i ± e =F , r ± e =F and r^fi^, depending on the forms of K', resulting in the low symmetry 
breaking scale of ~ TeV. Specifically, in order for the TeV symmetry breaking scale to occur, 
K' have to be in certain forms that will suppress the rare decays Kl — > /i ± e =F . In fact Ref. 
@ pointed that there are only 4 possible (approximate) forms for K' that are consistent 
with the TeV SSB scale: 



K' 3 




cos (3 sin (3 
1 

— sin (3 cos (3 _ 
/ cos 5 sin 5 

— sin 5 cos 5 

V i o o 



Note that Eq. ( p6|) is a special case of these forms, namely, 

K' 2 ((3 = 0)=K' 3 ( 1 = 0). 



(98) 



(99) 



Eq. ( p5|) amounts to suggesting that non-diagonal K' (required for low symmetry break- 
ing) and non-diagonal (required to obtain maximal — > v t l oscillations) may have a 
common origin. For example, if we make the ansatz that A3, A4 takes the approximate form 



(100) 



with the other 2 Yukawas {Ai,A2} approximately diagonal, then this will simultaneously 
lead to Eq. fl94]) and (p6|), hence the low symmetry breaking scale (~ few TeV) and near 
maximal — > i> t l oscillations. 
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B. Permutation of up- type quark masses and three generations mixing 

In this subsection we will include the first generation neutrino into the picture under the 
assumption that the intergenerational mixing between the first and the 2-3 sector is small. 
In such a scenario will take the approximate form 
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m L 



(-Ik \ 

g g r»22m33 

g "i22r»33 g 



101} 



'"R 

which leads to a pair of almost degenerate eigen masses m 2 , m 3 in the 2-3 sector that is 
separated from eigen mass 

mi = -^ = - meM ^ 2 eV (102) 

Ml 2m u 

with a distinct gap. The near maximal oscillation solution for atmospheric neutrino anomaly 
corresponds to a near degenerate pair of eigen masses m 2 ,m 3 . It is obvious that to include 
the first generation neutrino to solve the solar neutrino anomaly (which involves a much 
smaller 5m 2 ~ 10 -3 eV 2 required by all oscillation solutions) the first generation neutrino 
mass will have to be very nearly degenerate with the other two neutrino masses, which is 
quite unnatural. Thus we conclude that the gauge loop mechanism seems to explain the 
LSND data more readily than the solar neutrino anomaly. The gap 5mf 3 « 5m\ 2 can then 

be identified with the LSND measurement, 0.2 eV 2 ~ Sml SND ~ 3 eV 2 Q 

In the case of two generations discussed in the previous subsection, the form of the matrix 
K' = K' 2 {(3 = 0) = K' 3 (j = 0) in Eq. fl96|) leads to the scale of m 23 that are incompatible to 
Eq. fl92|). However, this result corresponds to only one specific form of K' in Eq. fl98l). In 
general, the other forms of K' in Eq. (|98|) could also lead to maximal neutrino oscillations 
that are consistent with TeV scale SSB, assuming that Yr = K' holds. Essentially there are 
only 4 special forms of K' that are of our interest. The case of K' = K' 2 (/3 = 0) = ^3(7 = 0) 
has been shown to be incompatible with Eq. (0). We will investigate the other three forms 
of K' in turns, namely (a) K' 4 (6 = 0) = K' 2 ((3 = f), (b) K[(a = 0) = K' 3 (j = |) and 
{c)K',{a = l) = K' A {5 = l). 

To find out how these three different forms of K' can also lead to maximal mixing in 
the 2-3 sector, we will use the constraint that the lop-sided form of Mr, and thus maximal 
v^l v t l oscillations, is preserved when takes on different form other than that of Eq. 
M as we replace Y ] R -> Y$. [Y$ are the other forms of K' as catalogued in Eq. fl98|) 1. 

Referring to Mr in Eq. (p8[) , since we know that, with the ansatz of Eq. (|88|) and the 
choice of basis, the matrices Mi, V, U', Me, Yl are diagonal, it then follows that the form of 
Mr goes like 



Mr ~ U^M U Y ] R ~ x . (103) 




W+Note that the solar neutrino problem can be solved in the current scenario if we mirror sym- 
metrize the model as we did in section VI so that maximal v e — * u' e oscillations result. In particular 
if we want to explore the possibility that the atmospheric neutrino anomaly is solved via near max- 
imal Un — > v T oscillations then this is actually compatible with mirror symmetry, since we just need 
to be in the parameter region where the oscillation length for — > v' oscillations is much greater 
than the diameter of the earth for atmospheric neutrino energies. 
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When is replaced by Y R via some transformation T, 

Yj -> f£ = TFjcj, (104) 
the lop-sided form of M R should be preserved, 

x 

U ] M U Y ] R -> U' ] M U TY ] R ~ | x ] . (105) 

0x0 



Since Y^ 



/l 0" 

1 [see Eq. (pi])], this mean that f/'f ought to conform to the condition 
\0 10. 



that 

U' ] M U T = I x I . (106) 




The matrix T relates the 'original' form of Y R [as in Eq. (|3)] to Yg via Eq. ( TO . Once 
we know T then we could work out from Eq. (|106|) . With Y R — > Y R and ?7^(= I) — > U'\ 
the net effect is that the diagonal up-type quark mass matrix M u is replaced by 

M u -> U^M U T, (107) 

thus changing the up-type quark masses dependence of the eigen masses m«. 

(a) Y£ = KM = 0) = K' 2 (0 = 

(0 1 0\ 

Let's look at the form of Y R = K' 4 (5 = 0) = K' 2 ((3 = f ) = 1 . The corresponding 

\i o o/ 

/0 1\ 

T matrix is T = 1 . The matrix U'\ by Eq. ( M , is C/'t = 

M u = diag {m u ,m c ,m t } -> U n M u T = diag {m t , m c , m u }. (108) 

As a result, the up-type quark masses as appear in m 2j 3 and mi in Eq. (|91~1) and (|102|) will 
be permuted by Eq. ( |108| ) which lead to the eigen masses for the light neutrinos 

fm T \ rfS 2 M E , fm e \ .„ . 

V m c / 2 \m ( / 

We see that the upper limit of the scale of m\ 3 are of the order 

m 2, 3 /eV 2 ~ 200?y 4 , (110) 

which permits a large range of parameter in rj so that m\ 3 ^> 5m 2 atrn for self consistency. 
The mass gap 
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m 9 



mi- 
me 



111) 



also permits a large range of parameter space in 77 to accommodate Smj^^. This scheme is 
easily realised by imposing the ansatz that the Yukawas are of the forms 



(112) 





fx 





0\ 




fo 





x\ 




/o 


X 


0\ 


A 2 = 





X 





, Ai = 





X 





, A 3i4 — 








X 




VO 





x/ 




V x 





0; 




U 





0/ 



(b) Y% = K[(a = 0) = ^(7 = f): 



Next, we look at the form 








[a = 


0) = ^(7 = 


7T 
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) = 


and T assume the forms 




















/o 


1 
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0\ 






1 





!) 


■ T -(i 












\° 












1/ 



1 

1 I . In this case, C/ /f 
1 



This results in permutation of up-type quark masses {m tt) m c ,m f } 
mass square difference 



1 2 2 1 

l m 2 — m ll 



ml 



(113) 

{m c ,m u ,m t }. The 
(114) 



is not compatible with the experimental values of Smj^^. In addition, the scale of m\ z is 
also too tiny to accommodate 5m 2 atm : 




A 



3,4 



m 2 2 3 < 4 x 10" 5 eV 2 . 
In passing, we note that the Yukawas of the forms 

fx \ 
A 2 = x , Ai 
V x) 

will lead to the above unsatisfactory scheme, 
(c) Yj = K' x {a = l) = K{5=1) : 

Finally, we investigate the case Y$ = K[(a — |) = K' A {8 
[/'' and T assume the forms 



(115) 
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(116) 
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In this case, 
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vo 
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0/ 



(117) 
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The up-quark mass permutation is {m u ,m c ,m t } — > {m t ,m u ,m c }. The upper bound of the 
scale of m\ 3 is 

m^/eV 2 £ 0.7r? 4 , (118) 

which means that it is possible to accommodate bmr atm . The mass squared difference for 
LSNDis 

Smlan = \ml - ml\ ~ m\ = SV< {^f Z 10" 2 (^eV 2 . (119) 

In this case we see that its parameter space can still accommodate the LSND and atmospheric 
neutrino anomalies [although the regime of parameter space is more restricted compared to 
the case (a)]. This scheme can be implemented if the Yukawas take the forms 

fx 0\ fO x\ /0x0\ 

A 2 = x , Ai = x , A 3i4 = x . (120) 

yoox; \xoo/ \x o o / 

In short, we see that out of the 4 different forms of K' in Eq. (|98|), the schemes (a) 
and (c) stand out to be most promising to provide viable solutions to both LSND and 
(near maximal oscillations) atmospheric neutrino anomaly with mass scale ~ few TeV. Note 
the schemes we have discussed in this section cannot accommodate solar neutrino anomaly 
(while keeping both LSND and atmospheric neutrino solutions) because the right-handed 
neutrinos are decoupled from the left-handed ones at TeV scale. 

In the parameter range where scalar interactions dominate over the gauge interactions 
things are less constrained and there are more possibilities. It is possible to implement the 
lop-sided Mr scheme in the scalar sector case in much the same way as in the previous 
section. Other possibilities where is off diagonal is also possible in the scalar case which 
can also lead to schemes compatible with data. However these schemes seem less elegant 
because of the larger degree of arbitrariness in scalar interactions. 



VIII. CONCLUSION 

The similarity of the quarks and leptons suggests that quarks and leptons might be 
connected by some spontaneously broken symmetry. However such a situation will lead 
to a gauge hierarchy problem unless the symmetry breaking scale is less that a few TeV. 
There are only two known ways that quark-lepton unification can occur at the TeV scale. 
First, quarks and leptons can be connected by a spontaneously broken discrete symmetry 
. While this an interesting possibility, it is difficult to naturally explain the lightness of the 
neutrinos in these schemes. The second possibility is a modification of the Pati-Salam model 
n called the alternative 422 model [§,0. It turns out that the neutrinos in the alternative 
422 model are naturally light because they are massless at the tree level and their masses 
are radiatively generated. The model also predicts novel B° physics. 

The possibility that the model can provide the interactions to generate the right neutrino 
mass and mixing patterns which might explain the atmospheric, solar and/or LSND neutrino 
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anomalies has been studied in detail in this paper. We have shown that the model cannot 
accommodate simultaneously all three of the anomalies unless it is extended in some way. 
However the minimal model can quite naturally accommodate the atmospheric and LSND 
anomalies. We have also pointed out that the solar neutrino problem could be most naturally 
explained if the model was extended with a mirror sector. 
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Appendix 

In this Appendix we will obtain the precise form of the would-be Goldstone bosons Gf R 

and physical Higgs fields Hf 2 in terms of the charged weak eigenstate fields \l r ■• an d 
0~5> - 5 [i. e. the matrix U g in the Eq. In addition we will also obtain the interesting 

result that 



lim V m 



lim V m 

"1^0 



0. 



(121) 



The Goldstone bosons associated with each spontaneously broken symmetry are given 



by 



G = $ T T a A a 



;i22) 



where $ are the Higgs scalar of the theory, T a the generators of the broken symmetries 
and A a the vacua. G\ is the would-be Goldstone bosons that are eaten by W L as their 
longitudinal polarisation. The associated (charged) generators of the broken SU(2)l are rf. 
Generally, 



Gf 



Xirt{XL) + Xl ^l(Xl) + XbTHxr) + Xr ^(Xr) + 

1 

~N~l 



c\T+, 



(123) 



For definiteness let we focus on the negatively charged fields, and let us work in the limit 
Mi 0: 



G] 



1 



wlXl 



1 _1 

' 2 ' 2 



(124) 



Likewise, Gf are associated with the SSB of charged sector in SU(2) R : 



G R 



XiTt(XL)+XLTt{xl)+XRTt{XR)+XRTt{X R ) + 

1 



<frU<P) + m i r^<P c ) 



(125) 



from which we obtain 



G 



R 



tt~> 



w 



R 



w RXr 



u 2 { 



1 1 

)2 ' 2 



r 



(126) 



The states orthogonal to Gf R will be the physical Higgs H^ 2 which in general not uniquely 
fixed because there are 2 directions orthogonal to Gf; R . 

A particularly important term in the Higgs potential is the term responsible for mixing 

4,-i 4,-i 

Xl 2 with Xr 2 which effects the generation of Dirac masses of neutrinos in the x-loop, 
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V 3 = Mxl<f>T 2X R + M'x\<\> c t 2X r + H.c. 



(127) 



We will show that there exists an approximate global symmetry U {l)x of this Higgs potential 
in the limit u\ — > 0, u 2 7^ (or u 2 — > 0, U\ ^ 0) that will allow us to identify the physical 
Higgs states (in that limit). The global symmetry U(l)x can be defined by the generator 



X = Y' + I aR - I 



<3L; 



where 



y'<xr) = -\> y'(xl) = \, Y'(<p) 



1. 



128) 



(129) 



Furthermore for the global symmetry to be useful we require it to be unbroken by the 
vacuum, that is, 



X( XL )=X( XR )=X^) = 



(130) 



which is indeed the case given our choice of X (and the limit U\ — *■ 0). 

Referring to V3 in Eq. (127), notice that the M' term is not a symmetry under U(l)x- 
However, in the limit u\ — > then M' must also be zero for self consistency. The reason 
is that a non zero M' in V3 will induce a linear term in the (21) component of <fi (i.e. the 
place where u\ would sit) when Xl,r develop VEVs. Because the potential is linear in 
u\ (for small Ui) a non-zero VEV for u\ must arise which is obviously not self consistent 
with our assumption that u\ = 0. This shows that U(l)x becomes an unbroken symmetry 
because the U(l)x asymmetry M' term vanishes in the limit u% — > 0. (We can also draw 
similar conclusion in the limit u 2 —>■ 0, U\ 7^ 0.) Furthermore it allows us to uniquely specify 
the physical Higgs fields if 1, if 2 since now we have two requirements. First they must be 
orthogonal to Gl,r and second they must be composed of components with the same U(l)x 
charge. (Note that ifi has X— charge -1, H 2 has X— charge +1.) These considerations lead 
to the identification of the charged physical Higgs fields: 



Hi 



Ho 



1 



1 + (u 2 /w R y 
1 



"2 



W R 



-Xr 



u 2 ( 



1 ll 

)2 ' 2 



u 2 \fi + (u 2 /w L y 



W L 



-Xl 



- U 2 (, 



;i3i) 



Writing Eq. ( 124 ), (|126|) , ( |131| ) in matrix form, we obtain Eq. (|43D. Note that the matrix 

because the rows (and columns) are all orthogonal. The existence 



U g is unitary, U~ 
of the approximate global symmetry X in the limit U\ 
means that the mass mixing term 



0,u 2 ± (or u 2 -> 0,«i ^ 0) 



V 3{4> ) oc uiu 2 . (132) 

The effect of this is that we can treat the mass mixing terms V^u) = —{Mu\ + 
4 — — 4 — —* 

M U2)Xl 2 Xr 2 as a small perturbation when u\ (u 2 ) is switched on from zero. 



30 



REFERENCES 



[1] J. Pati and A. Salam, Phys. Rev. D 10 (1974) 275. 

[2] R. Foot and H. Lew, Phys. Rev. D 41 (1990) 3502; Nuovo Cimento A104, 167 (1991); 

R. Foot, H. Lew and R.R. Volkas, Phys. Rev. D 44 (1991) 1531. 
[3] G. Valencia and S. Willenbrock, Phys. Rev. D50, 6843 (1994); A. V. Kuznetsov and N. 

V. Mikheev, Phys. Lett. B 329, 295 (1994); R. R. Volkas, Phys. Rev. D53, 2681 (1996); 

A. D. Smirnov, Phys. Lett. B431, 119 (1998). 
[4] H. Harari and Y. Nir, Nucl. Phys. B292, 251 (1987). 
[5] R. Foot, Phys. Lett. B 420 (1998) 333. 

[6] T. Haines et al., Phys. Rev. Lett. 57 (1986) 1986; Kamiokande Collaboration, K.S. 
Hirata et al, Phys. Lett. B 205 (1998) 416; ibid. B 280 (1992) 146; ibid. B 335 (1994) 
237; 1MB Collaboration, D. Casper et al., Phys. Rev. Lett. 66 (1989) 2561; R. Becker- 
Szendy et al., Phys. Rev. D 46 (1989) 3720; Nusex Collaboration, M. Aglietta et al., 
Europhys. Lett. 8 (1989) 611; Frejus Collaboration, Ch. Berger et al., Phys. Lett. B 227 
(1989) 489; ibid. B 245 (1990) 305; K. Daum et al, Z. Phys. C 66 (1995) 417; Soudan 2 
Collaboration, W.W.M., Allison et. al., Phys. Lett. B 391 (1997) 491; SuperKamiokade 
Collaboration, Y. Fukuda et al., Phys. Lett. B 433 (1998) 9; Phys. Lett. B 436 (1998) 
33; Phys. Rev. Lett. 81 (1998) 1562; MACRO Collaboration, M. Ambrosio et al, Phys. 
Lett. B 434 (1998) 451. 

[7] Homestake Collaboration, B. T. Cleveland et al., Astrophys. J. 496 (1998) 505; 
Kamiokande Collaboration, Y. Fukuda et al., Phys. Rev. Lett. 77 (1996) 1683; S. 
Fukuda, |hep-ph / 103032| ; SAGE Collaboration, J. N. Abdurashitov et al., Phys. Rev. 
Lett. 83 (1999) 4686; GALLEX Collaboration, W. Hampel et al., Phys. Lett. B 447 
(1999) 127; GNO Collaboration, M. ALtmann et al., Phys. Lett. B490, 16 (2000). 

[8] LSND Collaboration, W. Athanassapoulos et al., Phys. Rev. Lett. 81 (1998) 1774; Phys. 
Rev. C 58 (1998) 2489. 

[9] R. Foot and G. Filewood, Phys. Rev. D 60 (1999) 115002. 
[10] R. Foot, R.R. Volkas, O. Yasuda, Phys. Rev. D 58 (1998) 013006; P. Lipari, M. Lusignoli, 

Phys. Rev. D 58 (1998) 073005. 
[11] S. Fukuda et al, Phys. Rev. Lett. 85 (2000) 3999. 
[12] R. Foot, Phys. Lett. B 496 (2000) 169. 

[13] G. Conforto et al, Astropart. Phys. 5 (1996) 147 ; Phys. Lett. B 427 (1998) 314. 
[14] A. H. Guth, L. Randall and M. Serna, JHEP 9908 (1999) 018 . 
[15] R. Foot and R. R. Volkas, |hep-ph/95105T2l and also Ref. 
[16] R. Foot, Phys. Lett. B 483 (2000) 151. 

[17] R. Crocker, R. Foot and R. R. Volkas, Phys. Lett. B 465 (1999) 203. 
[18] R. Foot and R. R. Volkas, Phys. Rev. D 55 (1997) 5147. 

[19] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B 420 (1998) 397; Palo Verde 
Collaboration, F. Boehm et al, Phys. Rev. Lett. 84 (2000) 3764; Phys. Rev. D 62 (2000) 
072002. 

[20] Super-Kamiokande Collaboration, S. Fukuda, [hep-ex/0103032j . 
[21] Super-Kamiokande Collaboration, S. Fukuda, [hep-ex/0103033j . 
[22] SNO Collaboration, |nucI-ex/0l06TJT5" . 



see the Borexino website, [http: / /almime.mi.infn.it/ . 



31 



[24] see the KamLAND website jhttp: / /www. awa.tohoku.ac.jp/html/KamLAND/ 



[25] J. Bunn, R. Foot and R. R. Volkas, Phys. Lett. B 413 (1997) 109; R. Foot, R. R. Volkas 

and 0. Yasuda, Phys. Rev. D (1998) 1345. 
[26] |http:/ /www. neutrino.lanl.gov/BooNE/ 



[27] G. Raffelt and D. Seckel, Phys. Rev. Lett. 60 (1988) 1793. 

[28] R. Barbieri and R. N. Mohapatra, Phys. Rev. D 39 (1989) 1229. 

[29] P. Langacker and S. U. Sankar, Phys. Rev. D 40 (1989) 1569. 

[30] V. Berezinsky, in The Ninth International Workshop on Neutrino Telescopes, Venice, 

^stro-ph/010731161 (2001). 
[31] M. S. Tuner, Phys. Rev. Lett. 60 (1988) 1797. 
[32] K. Kanaya, Prog. Theor. Phys. 64 (1980) 2287. 

[33] T. L. Yoon and R. Foot, Phys. Lett. B 491 (2000) 291; B. C. Allanach, Phys. Lett. B 
450 (1999) 182. 

[34] R. Foot, H. Lew and R. R. Volkas, Phys. Lett. B 272 (1991) 67; R. Foot, H. Lew and 
R. R. Volkas, Mod. Phys. Lett. A 7 (1992) 2567; R. Foot, Mod. Phys. Lett. A 9 (1994) 
169; R. Foot, R. R. Volkas, Phys. Rev. D 52 (1995) 6595. 

[35] R. Foot, |astro-ph/0102294] . 



32 



